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ABSTRACT 

We present optical photometry of the X-ray transient XTE J2 123-058, obtained in 
July-October 1998. The light curves are strongly modulated on the 5.95hrs orbital 
period, and exhibit dramatic changes in amplitude and form during the decline. We 
used synthetic models which include the effect of partial eclipses and X-ray heating 
effects, to estimate the system parameters, and we constrain the binary inclination 
to be i=73° ± 4. The model is successful in reproducing the light curves at different 
stages of the decay by requiring the accretion disc to become smaller and thinner by 
30% as the system fades by 1.7 mags in the optical. From Aug 26 the system reaches 
quiescence with a mean magnitude of R=21.7±0.1 and our data are consistent with 
the optical variability being dominated by the companion's ellipsoidal modulation. 
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1 INTRODUCTION 

Soft X-ray transients (SXTs) are a subclass of low-mass X- 
ray binaries (LMXBs) that are characterized by episodic X- 
ray outbursts (usually lasting for several months) , when the 
X-ray luminosities can increase by as much as a factor of 
10^ (van Paradijs & McClintock, 1995). The observed opti- 
cal flux is generated by X-ray re-processing in the accretion 
disc and the companion star. These outbursts recur on a 
time scale of decades, but in the interim the SXTs are in 
a state of quiescence and the optical emission is dominated 
by the radiation of the faint companion star. This offers the 
best opportunity to analyze the properties of this star and 
obtain dynamical information which eventually enables us to 
constrain the nature of the compact object. There are cur- 
rently 12 SXTs with identified optical counterparts, with 8 
dynamical black-holes and 3 confirmed neutron stars: CenX- 
4, Aql X-1 and 1608-522 (van Paradijs & McClintock, 1995). 
In addition, there are a few neutron star binaries exhibit- 
ing X-ray on and off states ( EXO0748-676, 4U 2129+47 
and SAX J1748.9-2021 ), although they are not classified 



as SXTs because they do not show the classic fast rise and 
slow exponential/linear decay. 

The X-ray transient J2123-058 was discovered on 29 
June 1998 by the Rossi X-Ray Time Explorer (RXTE) 
(Levine et al. 1989) reaching a peak X-ray flux of 100 
mCrabs (2-12 keV). Its high Galactic latitude (b=-36°.2) 
is unusual among transients, an indication that J2123-058 
might be a member of the galactic halo population. The 
optical counterpart was promptly identified with a vari- 
able star of R=17.2 (Tomsick et al. 1998a), which was only 
marginally visible on a digitized U.K. Schmidt plate, sug- 
gesting a preoutburst magnitude R>20 (Zurita et al. 1998). 
Spectra obtained early in the outburst showed strong high- 
excitation lines of Hell A4686, Ciii/Niii A4640 and weak 
Balmer emission embedded in broad absorptions (Tomsick 
et al. 1998a, Hynes et al. 1999). These features are fre- 
quently observed in SXTs during outburst (e.g. Callanan 
et al 1995) and persistent LMXBs (e.g. Augusteijn et al. 
1998). Type-I (thermonuclear) bursts have been detected 
both in X-rays (Takeshima and Strohmayer 1998) and op- 
tical (Tomsick et al. 1998b), a signature of a neutron star 
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in the binary. The outburst light curve exhibited regular 0.7 
mag deep triangular-shaped minima repeating every 6hrs 
(Casares et al. 1998, Tomsick et al. 1998b), a strong indi- 
cation of high inclination. This provided the first evidence 
for the system orbital period (P= 5.957±0.003) which was 
later confirmed by a radial velocity study of the Hell A4686 
emission line (Hynes et al. 1998). In addition, Ilovaisky and 
Chevalier (1998) reported the presence of a 0.3 mag mod- 
ulation with a period of 7.2 days, probably caused by the 
prcccssirig disc. Since 26 Aug the system had settled down 
to its quiescent state at R~21.7 (Zurita & Casares 1998). 

This paper presents the results of a comprehensive set 
of observations that have led to detailed optical light curves 
from outburst through the decay into quiescence. Our spec- 
troscopy will be the subject of a second paper (Hynes et al. 
2000). 
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^lACSO-SOcm Telescope in the Observatorio del Teide (Tenerife). 
^ 1 m Optical Ground Station in the Observatorio del Teide 

(Tcncrifc). 

^ Mount Canopus Im telescope at Tasmania. 

^ 0.5m telescope at Crimea 

^ 1.25m reflector of SAI Crimean Station. 

^ 1.2m telescope at the Kryonerion Astronomical station of the 
National Observatory of Athens at Kryonerio of Korinthia. 
^ Im Jacobus Kapteyn Telescope in the Roque de los Muchachos 
Observatory (La Palma). 



2 OBSERVATIONS AND DATA REDUCTION 

We observed J2123-058 during the period 1998 July- 
September with the IAC80 and Im Optical Ground Station 
(OGS) in the Observatorio del Teide (Tenerife), the 1.25m 
SAI Station and 0.5m and 0.38m telescopes at Crimea, the 
Mount Canopus Im telescope at Tasmania and the 1.2m 
telescope at the Kryonerion Astronomical station of the Na- 
tional Observatory of Athens at Kryonerio of Korinthia. On 
the night of 21 of June of 1999, we obtained VRI images, us- 
ing the Im Jacobus Kapteyn Telescope in the Roque de los 
Muchachos Observatory (La Palma) . The integration times 
ranged from 1 to 40 min, depending on the telescope, at- 
mospheric conditions and the brightness of the target. The 
observing log is presented in Table 1. All images were cor- 
rected for bias and flat-fielded in the standard way. 

We applied aperture photometry to our object and sev- 
eral nearby comparison stars within the field of view, using 
IRAF. We selected four comparison stars which were checked 
for variability during each night, and over the entire data 
set. We calibrated the data was using 9 standard stars from 
several fields (Landolt 1992), from which we constructed a 
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colour dependent calibration. The Crimean data in R were 
calibrated relative to the nearby comparison star from the 
USNO catalogue which was also in to the previous calibrated 
data. White light photometry was also obtained in Crimea 
using a 0.5m telescope equipped with a TV detector. The TV 
camera operates in the range AA3500-8000 with maximum 
sensitivity at A4000 and effective wavelength at Ae//4851. 
Magnitudes were measured relative to a neaxby star from 
the USNO catalogue which was calibrated at B and V wave- 
lengths. 'Equivalent BV magnitudes', were calculated by in- 
terpolating the fluxes in B and V to find the flux at the 
effective wavelength of the TV detector. For further details 
see Pavlenko, Prokofieva & Dolgushin (1989). 



3 LIGHT CURVES 

3.1 Long term light curve 

Figure 1 compares the overall light curve of J2123-058 
in optical (VR bands) and X-rays (2-12 keV). The long- 
term behaviour in X-rays shows a classical FRED (fast-rise 
exponential-decay) morphology, with characteristic e-folding 
times of 2.4 d (rise) and 19 d (decay). These time scales 
coincide with the mean values of the distributions of SXTs 
(Chen, Shrader & Livio 1997). Twenty days after the peak of 
the outburst, the X-ray intensity reaches a secondary maxi- 
mum (of less than half of the peak intensity) . The secondary 
maximum is also suggested by our optical data. 

We identify 3 different stages in the optical light curve: 
the outburst plateau (until ~ 10 Aug), the decay phase (~ 
10-26 Aug) and quiescence (from ~ 26 Aug). In the plateau 
phase the object brightness decays at a moderately slow rate 
of ~ 0.03 mag day"^ although a modulation is clearly visi- 
ble in the nightly mean magnitudes. The time scale of this 
variability is consistent with the 7-d modulation attributed 
to disc precession by Ilovaisky & Chevalier (1998). From ~ 
10 Aug the optical light curve began an abrupt fall at a rate 
of ~ 0.2 mag day^^ before reaching quiescence on 26 Aug 
at R=21.7 (sec Figure 1). 

Taking V(pcak)~17.2 and V(quiesccnt)~22.9, we esti- 
mate a total outburst amplitude of 5.7 mags. Using the em- 
pirical relation AV=14.36-7.6- lgPori.(hr) (Shahbaz & Ku- 
ulkers 1998), we would expect a total amplitude AV~8.4, 
2.7 mags larger that what is observed. This difference can be 
explained by 3 effects: (1) the outburst brightness of the disc 
being reduced by a factor cos(i), since the disc would be fore- 
shorted by the high binary inclination angle, (2) we are also 
assuming that the quiescent flux is completely dominated by 
the companion star with no veiling from the accretion disc 
and (3) if the secondary star is sufficiently evolved for it to 
be degenerate it would appear intrinsically fainter. We be- 
lieve the discrepancy is due to a combination of these effects 
and can only be resolved by obtaining optical spectroscopy 
in quiescence. 

We see evidence for optical bursts both during outburst 
(29 July with an amplitude of 0.3 mag) and at the onset of 
quiescence (27 Aug and 2 Sept with amplitudes larger than 
1 mag). The optical bursts observed during the ourburst 
decay is most probably caused by re-processing of X-ray 
bursts. However, the origin of the optical bursts observed 
during the onset of quiesence is more puzzling, given that 



the source had almost reached its quiescent (low luminosity) 
X-ray state. Note that no simultaneous X-ray observations 
exist, which might shed light on the origin of these optical 
flashes observed during quiescence. 

Our V and 'BV (Figure 1), show that the BV magni- 
tudes drop faster than R. Moreover the amplitudes in BV 
and R increase during the decay. Figure 2 presents colour 
information of J2123-058 as a function of time. Although 
we cannot directly compare the 'BV'-R colour with V-R it 
is clear from the plot that the system redden as the out- 
burst decays and the secondary's contribution increases. In 
quiescence, we obtained an upper limit to the quiescent V 
magnitude through 2x2400s images of ,12123-058 with the 
IAC80 telescope. We also have a marginal detection in V 
on the night of 21 Jun 1999 using the JKT telescope. Our 
colours are consistent with the spectral type of a late-K 
main sequence star. 

3.2 The orbital light curve 

Representative light curves of J2123-058 in different stages 
of the outburst cycle are presented in Figure 3. The data 
have been folded on the updated ephemeris given by Zurita 
& Casares (1998): HJD 2451042.639(5) + 0.24821(3) E. Note 
the dramatic changes in amplitude and morphology of the 
light curves as the outburst decays. 

The July light curves are flat topped with broad trian- 
gular minima. They have a full amplitude of 0.7 mag and 
are reminiscent of the 5.1 hr period eclipsing transient EXO 
0748-676, although its long term behaviour is characterized 
by high and low X-ray/optical states rather than transient 
outbursts. We also find similarities with the eclipsing LMXB 
source 2A1822-371 whose complex optical light curve has 
been successfully modeled by X-ray heating and partial oc- 
cultation of a thick non-axisymmetric accretion disc and a 
faint companion (Mason et al. 1980). Asymmetries in the 
eclipse minima, are also seen in J2123-058 on individual 
nights and also in the curves by Tomsick et al. 1998. 

On the other hand, the light curve of Aug 16 is almost 
sinusoidal and shows a peak-to-peak amplitude of ~ 1.4 mag. 
It resembles the outburst light curve of the 5.2 hr neutron 
star system 4U 2129-1-47, where the large amplitude has been 
attributed to X-ray reprocessed radiation from the heated 
face of the optical star (Thorstensen et al.l979). Note, how- 
ever, that at the onset of the fast optical decay (~10 Aug) 
the X-ray flux had already dropped by a factor > 10 with 
respect to the outburst peak. Note also the presence of two 
narrow 0.2 mag dips at phases and 0.5 which suggest the 
presence of grazing eclipses. 

Finally, the quiescent light curve displays a characteris- 
tic ellipsoidal modulation from the secondary star: a double- 
humped variation with a full amplitude of ~0.4 mag. This 
light curve was produced by phase binning the entire quies- 
cent data (from Aug 26 on) after detrending the night-to- 
night variability using a linear fit. 



4 DISTANCE ESTIMATE 

In the context of King & Ritter's (1998), the X-ray expo- 
nential decay seen in J2123-058 indicates that irradiation 
is strong enough to ionize the entire accretion disc. Also, a 
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secondary maximum is expected one irradiated-state viscous 
time after the onset of the outburst and it can be used to 
calibrate the peak X-ray luminosity and hence the distance 
to the source Dkpc through 

Okpc = 4.3 X W-h'JV^^f^^F^-'^\,'/^ 

where Fp is the peak X-ray flux, ta the time of the secondary 
maximum after the peak of the outburst in days, Td the e- 
folding time of the decay in days, rj the radiation efficiency 
parameter and / the ratio of the disc mass at the start of the 
outburst to the maximum possible mass (Shalibaa, Charles 
& King 1998). In our case, r£j=19 d, tg ~20 d and Fp can 
be estimated from the XTE count rate (6 counts s^^ in the 
energy range 2-12 keV) which corresponds to 1.84 x 10~® 
erg cm~^ . Assuming 77=0.15 and /=0.5 we find Dkpc = 
5.7. 

Alternatively, we can estimate the distance to the source 
by comparing the quiescent magnitude with the absolute 
magnitude of a main sequence star which fits within the 
Roche lobe of a 6hr period orbit. Combining Paczynski's 
(1971) expression for the averaged radius of a Roche lobe 
with Kepler's Third Law we get the well-known relationship 
between the secondary's mean density and the orbital pe- 
riod: p = llO/Phj. (g cm~^). Substituting the orbital period 
of J2123-058 we obtain p=3.1 g cm~^ which corresponds 
to a K7V secondary star with mass ~0.6 M© and absolute 
magnitude Mr ~ 7. The dereddened quiescent magnitude 
is R=21.4 (using Av=0.37±0.15 as derived from the NaDl 
line; see Hynes et al. 1998) which yields Dkpc = 7.7. Strictly 
speaking, this is a lower limit to the distance as we are ne- 
glecting any contribution by the accretion disc to the qui- 
escent optical flux. However, note that the true distance is 
probably not too far off 8 kpc since our quiescence light curve 
does not show strong evidence for disc contamination (see 
section 5 and Figure 6). A conservative limit can therefore be 
provided by assuming a 50 percent contribution by the ac- 
cretion disc to the continuum light. Allowing for 50 percent 
disc contamination (as observed in J0422-I-32, a black-hole 
transient with comparable orbital period; see Casares et al. 
1995) we obtain R=22.2 for the companion star and hence 
-Dkpc = 10.8. Hereafter we will adopt Dkpc = 8 ± 3 which is 
consistent, at the lower end, with other estimates based on 
photospheric expansion models of the X-ray bursts (Homan 
et al 1999, Tomsick et al 1999). A spectral type determina- 
tion of the companion star is essential to refine this distance 
estimate. 



and the companion star are fixed to 0.95 and 0.40 respec- 
tively, following the results of de Jong et al (1996). The 
model parameters are the binary inclination (i), mass ratio 
(q=Mi/M2), the accretion disc radius (Rdisc) defined as a 
fraction of the distance to the irmer Lagrangian point (Rli), 
the flaring angle of the accretion disc (a) and the X-ray lu- 
minosity (Lx). 

In order to model the outburst light curve we have aver- 
aged our best quality light curves of the plateau phase (26- 
30 ,hily) in 29 phase bins. Using La; = 1.3x 10''^ erg s^^ (for 
-Dkpc = 8) and M2=0.6 M0 we performed a least-squares 
fit to the data. Our best fit solution gave a reduced chi 
squared xl of 1.36 for ?:=76.0±1.0 degrees, Rdisc=0.75i;;:;;^ 
Rli, a=7.6lQ'2 degrees and g=4.6tg"2. The uncertainties 
quoted are at the 99 percent confidence level and were ob- 
tained by grid searching the parameter of interest whilst op- 
timizing the other model parameters. We have also rescaled 
the Xv values so that the minimum is 1. The best model 
fit to the outburst data is shown in Figure 4 (solid lino). 
In order to examine the effects of changing M2, we fitted 
the outburst data with M2=0.1Mo (Figure 4 - dotted line). 
We find that the derived parameters are the same, within 
the errors: xl = 1-5, i=75.3't°;l degrees, Rdisc=0.83±0.05 
Rli, g=3.4±0.4, a—8.8'ti[l degrees. Regarding the decay 
stage, we have fitted the light curve of 16 Aug binned into 
20 phase bins, when the X-ray luminosity was an order of 
magnitude lower, La;=1.3xl0^^ erg s~^. We fitted the decay 
stage data with this X-ray luminosity and M2=0.6Mq. Our 
best fit solution gave a reduced Xi' of 1.51 for i=72.0±3.0 de- 
grees, Rdisc=0.56±0.06 Rli, a=5.7±0.5 degrees and g=4.2 
{q could not be constrained with this data set). The best 
model fit to the decay stage data is shown in Figure 5. Again 
99 per cent confidence levels arc quoted. 

In Figure 6, we show the quiescent light curve which 
exhibits the characteristic ellipsoidal modulation of the sec- 
ondary star: i.e. two equal maxima and two unequal minima. 
The miriimun at phase 0.5 is expected to be deeper than the 
minimum at phase 0.0 because gravity darkening is more 
pronounced near the inner Lagrangian point Li . This effect 
is important for high inclination systems (see eg. Avni and 
Bahcall, 1974). In the figure we also show model plots using 
1=73°, q=4.6 (values which are consistent with those derived 
from the decay and outburst data), and no X-ray heating. 
The solid and dotted lines show plots with zero and 50 per- 
cent disc contamination. The former model probably best 
describes the data. 



5 MODELLING 

In an attempt to interpret the different light curves and de- 
rive the system parameters, we have used a model based 
on the work by de Jong et al. (1996, also see references in- 
cluded). The model assumes a flared accretion disc of the 
form h (X r''^^ (with h and r the disc height and radius re- 
spectively) and a Roche lobe filling secondary and accounts 
for X-ray heating, shadowing effects and mutual eclipses of 
the disc and the secondary. The disc is assumed to radi- 
ate as a blackbody with a radial temperature distribution 
calculated according to Vrtilek et al. (1990). The intensity 
distribution on the secondary star is computed using Ku- 
rucz model atmospheres. The albedo of the accretion disc 



6 DISCUSSION 

XTE J2123-058 is a remarkable neutron star binary. Our 
optical light curve shows marked orbital modulations with 
dramatic variations as the outburst declines. These are 
very similar to the modulations observed in accretion disc 
corona (ADC) sources of comparable orbital periods (e.g. 4U 
2129+47 and 2A 1822-371; see eg. McClintock et al. 1980 
and Mason et al. 1982) although none of these are transient 
systems. 

The X-ray light curve displays the classic SXT proper- 
ties, namely a FRED morphology with typical e-folding rise 
and decay times and secondary maximum. Furthermore, the 
ratio of X-ray to optical luminosity = -Bo -1-2.5 log Fx(/uJy)] 
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is also in excellent agreement with the observed distribu- 
tion of LMXBs. Taking B=17.28 (Tomsick et al. 1998a) and 
Fx(2- 12fceF) ~ 100 mCrab (Levine, Swank & Smith 1998) 
at the outburst peak and assuming — 0.37 (Hynes et al 
1998) we obtain ^=21. 9, whereas the distribution peak of 
LMXBs gives 5=21.8 ± 1.0 (see van Paradijs & McClintock 
1995). This result implies that, despite its high inclination, 
the X-ray source in J2123-058 is not hidden by the accre- 
tion disc (i.e. a >90°-i). This is consistent with the values 
our model favors for the binary inclination (i=73°) and disc 
flaring angle (a=5°.7 - 7°. 6). 

The longterm evolution of the optical light curve can 
be compared to those of other SXTs (e.g. GRO J0422-I-32, 
A0620-00, N. Muscae 91). They show a slow linear decay 
followed by a steeper fall. We find that J2123-058 also re- 
produces this behaviour, although the total amplitude and 
time scales are a factor of ~ 2 shorter (see e.g. Callanan et 
al 1995). 

We have modeled our R-band light curves of J2123- 
058 at different stages of the outburst including the ob- 
scuration effects and X-ray heating of the secondary star 
accretion disc. This led us to constrain the system inclina- 
tion to i=73° ± 4. We find encouraging the excellent agree- 
ment between the inclination values obtained for the two 
independent light curves (July 26-30 and 16 Aug) during 
which Lx has dropped by one order of magnitude. The light 
curves at the plateau phase (July) are very similar to those 
of EXO 0748-676 and 2A 1822-371 with an extended de- 
pression of the luminosity from phase ~0.7 until the eclipse 
and a steeper rise to maximum (see Mason et al. 1980 and 
Schmidtke et al 1987). Our model fit indicates that the ac- 
cretion disc is the dominant source of light and the triangular 
shaped minima can be interpreted as eclipses of the accre- 
tion disc by the secondary star together with the changing 
aspect of the heated polar caps of the companion star. The 
dramatic changes observed in the light curves during decline, 
are triggered by large changes in the disc size and geome- 
try. Our fit to the decay data (16 Aug) demands a thin- 
ner and smaller accretion disc which implies a smaller frac- 
tion of the disc is X-ray heated. Conversely, the secondary 
star is more exposed to the X-ray radiation and therefore 
the total amplitude of the modulation increases by a fac- 
tor of 2 (to ~ 1.4 mag). The resulting light curve has a 
sinusoidal-like shape and is reminiscent of the LMXBs 4U 
2129+47(=V1727 Cygni) and HZ Her. Our model fits im- 
plies a change of ~30 percent in the disc size, as the system 
fades by 1.7 mags in the optical. The change in the radius 
of the disc size is what one expects. If angular momentum is 
transported outwards in the disc through viscous processes, 
then at outburst, since matter diffuses inwards, the angular 
momentum of that matter has to be transfered to the outer 
parts of the disc, and the radius of the disc is expected to 
expand. When the system is decaying, after the end of the 
mass transfer enhancement, the disc shrinks to its original 
radius (Livio & Verbunt 1988; Ichikawa & Osaki 1992). Ob- 
servations of U Gem, Z Cha, OY Car, HT Car show that 
the accretion discs are indeed larger in outburst than in 
quiescence (Smak 1984b; O'Donohuge 1986; Harrop-AUin & 
Warner 1996). Comparing our results with disc radius vari- 
ation in U Gem (Smak 1984a), we find approximately the 
same rate ofdecrease. 
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Figure 1. Left: Temporal evolution of J2123-058 plotted as 'X-ray magnitudes'(top) and R-band magnitudes averaged per day (bottom). 
Note that 1 Crab equals an asm count rate of 75 counts per sec. The X— ray data arc provided by the asm / rxte teams at mit and at the 
RXTB SOP at NASA's GSFC. Right: Temporal evolution plotted as R-band magnitudes (top) and V and 'BV magnitudes (bottom). 
The optical data points have been obtained with the following telcscopcs/sites: OGS (triangles), IAC80 (filled circles), Crimean telescopes 
(open squares), Kryoneri Astronomical Station (diagonal crosses) and Tasmania (open circles). Asterisks mark the magnitudes reported 
by Tomsick et al. (lAUC 6957). 
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Figure 2. (BV-R), (B-V) and (V-R) colours of J2123-058 as a 
function of time. Indices (B— V) and (V-R) were obtained from 
single simultaneous points. (BV-R) were calculated averaging 
magnitudes simultaneous within 3 min. Horizontal lines mark (V— 
R) colours for different spectral types. 



Figure 3. 

Optical light curves of J2123-058 taken during July-Sep 1998, 
covering outburst, decay and quiescence (bursts have been 
suppress in the averaged curves). They show spectacular 
variations, evolving from a strong triangular shape to sinusoidal 
and double-hump modulations. The light curves are folded on 
the 5.95 hr period and are shown twice for clarity 
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Figure 4. Outburst light curve (Jul 26-30) (points) and 
model fits (solid line) for La,=1.3xl037erg s"!, i=76°, q=4.6 
(M2=O.6M0), Rdisc=0.75 Rli and a=7°.6. The dotted line 
shows the best fit using M2=O.1M0 (see text). 
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Figure 5. Decay light curve (16 Aug) (points) and 
model fits(solid line) for La;=1.3xlO=^'3gj.g g-i^ i=72°, q=4.2 
(M2=O.6M0), Rdisc=0.6 Rli and a=5°.7 
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Figure 6. Quiescent light curve (26-27 Aug 98 and 24-26 Sep98) 
(points) and model plot (solid line) for i=73°, q=4.6 and no X- 
ray heating. The dashed line shows a model plot assuming 50 per 
cent contamination by the accretion disc. 



© 0000 RAS, MNRAS 000, 000-000 



